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Abstract 

We present the results of kinetic modelling of quantum effects in laser-beam interaction. In the developed numeri- 
cal model, electron-positron pair production by hard photons, hard photon emission and the electromagnetic fields 
generated by the created charged particles are taken into account. Interaction of a relativistic electron beam with a 
strong laser pulse is analyzed. It is shown that the quantum eff'ects can be important even for moderately intense laser 
pulses when the number of emitted photons by single electron is not large. Electron-positron pair plasma production 
in extremely-intense laser field via development of electromagnetic cascades is also studied. The simulation results 
confirm the prediction of strong laser field absorption in the self-generated electron-positron plasma. It is shown that 
the self-generated electron-positron plasma can be an efficient source of energetic gamma-quanta. 



Electromagnetic cascades are one of the basic phe- 
nomena of strong-field physics. As opposed to clas- 
sical physics, in quantum physics the decay of high- 
energy photon in strong electromagnetic fields is pos- 
sible. The photon can decay with creation of electron- 
positron (e~e'^) pair. If the energy of the electron and the 
positron is high enough, they can emit new photons that 
can, by-turn, decay with creation of new e~e'^ pairs, etc. 
This cascade process is also called 'electromagnetic (or 
QED) shower' . 

Due to its importance, the cascades recently are in- 
tensively studied. For example, electromagnetic show- 
ers that develop during the pass of high-energy particle 
(electron, positron or photon) through constant homo- 
geneous magnetic field are theoretically investigated in 
Ref. [1, 2]. In this case the energy for creation of sec- 
ondary particles is transferred from the initial energy of 
the seed particle. In such shower the mean energy of 
a particle decreases with the lapse of time. The dy- 
namics of electromagnetic showers in strong fields of 
pulsars also attracts much attention |3]. Showers in in- 
homogeneous, varying electromagnetic fields (fields of 
pulsars, for example) can be self-sustained. This means 
that secondary particles in such shower are accelerated 
by strong external electric field and the mean energy of 
a particle in such cascade is constant or increases. 

There are little experiments devoted to investigation 
of self- sustained cascades. However, the recent progress 
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in laser technologies will allow to fill up this gap. The 
laser intensity can be very high in near future so that 
quantum electrodynamic (QED) effects will be essential 
14 lil H B Ish • An electron (positron) moving in strong 
laser field can be accelerated up to very high energy and 
can produce electromagnetic cascade. If the laser inten- 
sity is high enough, the number of produced pairs can 
be so great that self-generated electromagnetic fields of 
e~e'^ plasma can strongly affect the further cascade dy- 
namics. Moreover, a large portion of the laser energy 
can be absorbed by self-generated e~e'^ plasma. This ef- 
fect is especially important to determine the limitations 
on the intensity of high power lasers |6,^. The estima- 
tions of intensity threshold for cascade development in 
circularly polarized standing electromagnetic wave are 
presented in Ref. L6]. However, in order to study the 
cascade development and to verify estimations the self- 
consistent numerical models is needed. 

We develop the self-consistent two-dimensional nu- 
merical model based on particle-in-cell (PIC) and 
Monte Carlo (MC) methods |9]. The model uses 
the probability rates of photon emission and electron- 
positron pair production calculated in the framework of 
QED theory 113, El El- It exploits the strong differ- 
ence between the photon energy of laser-plasma fields 
and the characteristic energy of hard photons emitted by 
relativistic charged particles in electromagnetic fields. 
Making of use standard PIC technique the hard pho- 
tos, electrons and positrons are modeled as quasiparti- 
cles while the laser and plasma field is calculated by nu- 
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Figure 1: Energy spectrum of the electron bunch after its collision 
with a linearly polarized laser pulse characterised by = 5 and ctx - 
32. The gamma-factor of the bunch electrons before the collision is 

ro = 2-io3. 



Figure 3: Energy spectrum of the electron bunch after its collision 
with a linearly polarized laser pulse characterised by = 100 and 
cTjc = 1.6. The gamma-factor of the bunch electrons before the colli- 
sion is 70 = 2 • 10^. 
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Figure 2: Energy spectrum of the electron bunch after its collision 
with a linearly polarized laser pulse characterised by ao = 20 and 
(Tx = 1.6. The gamma-factor of the bunch electrons before the colli- 
sion is 70 = 2 • 10^. 



merically solving of Maxwell's equations. The model 
has been benchmarked to the simulations performed by 
other MC codes 1 13]. The PIC part of the model is two- 
dimensional version of the model used in Ref . 

One of methods to probe QED in laboratory condi- 
tions is to study the interaction of a relativistic electron 
beam with an intense laser pulse |15, 16]. In the rest 
frame of the relativistic particle the laser field is very in- 
tense and can be close to critical. So at this field strength 
quantum eff'ects become important. In this case the en- 
ergy of a photon emitted by a beam electron undergo- 
ing oscillations inside the laser pulse can be close to the 
electron energy. 

We simulate the interaction between a relativistic 
electron bunch and a laser pulse by the means of our 
numerical model. The laser pulse is linearly polarized, 
has the Gaussian envelope Ey = = ao exp[-y^/cry - 



(x - t)^/crl], where E and B are the electric and mag- 
netic field envelopes, respectively, gq = \e\Eo/(mcoj), 
Eq is the maximal field strength of the laser pulses, 
oj = 2nc/A, A is the laser wavelength, m and e are the 
electron mass and charge, respectively, c is the speed 
of light, (Tx is the pulse duration, cTy is the pulse width. 
From here on coordinates are normalized to A and time 
is normalized to In/oj. The parameters of the simula- 
tion are the following. The wavelength /I = 0.8 yum, 
(Tx = 32, ao = 5, the laser spot size o-y is assumed to be 
much large then the bunch width. The initial gamma- 
factor of the bunch electrons is yo = 2000. For these pa- 
rameters;^ ^ 0.03 ^ 1 and the photon emission regime 
should be classical. The energy distribution of the elec- 
tron bunch after passing through laser pulse is shown 
in Fig. [U The mean gamma-factor of the electron after 
interaction is y = 1568. The task about the interaction 
between an electron and a laser pulse has exact theo- 
retical solution in the framework of Landau - Lifshitz 
representation of radiation reaction force 1 17]. This so- 
lution yields the following gamma-factor of the electron 
passing the laser pulse with Gaussian envelope: 
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(1) 



where ji = 2e^oj/(3mc^). For the given parameters this 
formula yields y = 1580 that is in good agreement with 
the simulation result. 

It is interesting to note that quantum nature of pho- 
ton emission can reveal itself even in classical regime 
(X ^ 1). It is seen from Fig. [T] that the energy spread 
of the electron bunch is large in contrast to the classical 
theory. The reason is that each bunch electron emits a 
small number of photons. The characteristic time be- 
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Figure 4: At the instant of time t ^ 24, spatial distribution of (a) 
electron and (b) photon normalized densities pe^y - ne^yliaQricr) and 
(c) laser intensity pi = (E^ + B^)/Eq averaged over laser wavelength 
during the collision of two laser pulses with = 1-5 • 10^. 
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Figure 5: At the instant of time t ^ 54, spatial distribution of (a) 
electron and (b) photon normalized densities pe,y = ne,y/(aoncr) and 
(c) laser intensity pi = {E^ + B^)/Eq averaged over laser wavelength 
during the collision of two laser pulses with «o = 1-5 • 10^. 
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Figure 6: At the instant of time t ^ 62, spatial distribution of (a) 
electron and (b) photon normalized densities pe,y = ne,y / (aoricr) and 
(c) laser intensity pi = {E^ + B^)/E^ averaged over laser wavelength 
during the collision of two laser pulses with = 1-5 • 10^. 



tween consecutive photon emissions (photon emission 
time) can be estimated as the ratio of the characteris- 
tic energy of emitted photon to the radiation power . In 
the classical limit the photon energy is hcj^m ~ aoy^hco 
and in the quantum limit hojem ~ ymc^ Uul . Making of 
use the expression for the radiation power in the classi- 
cal and quantum limits [HI] we can find the photon 
emission time in the classical and quantum limits 



(2) 



where ctf is so-called radiation formation length, which 
is equal to length of the electron trajectory path, over 
which the particle is deflected by angle 1 /y 111, 18.. 191. 
The radiation formation length for an electron collid- 
ing with ultra-high intense laser pulse can be estimated 
as follows: tf = l/(2nao). For the given parameters 
o'x/ictrad) - 6 that is not much larger than unity. There- 
fore the number of photons emitted by a bunch electron 
passing through the laser pulse is not large. Because 
of quantum nature of photon emission the spread in the 
number of the photons and energy of the photons emit- 
ted by bunch electrons can be significant. 

The electron energy distribution of the electron bunch 
after interaction with shorter and more intense laser 
pulse (cTx = 1.6, Go = 20) is shown in Fig. [2l The 
initial electron gamma-factor is 70 = 2000. For these 
parameters;^ ^0.1 and trad - l-l- As (Txltrad - 1 a 
bunch electron emits approximately one photon in the 
average. As the electron distribution function peaks at 
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Figure 7: At the instant of time t ^ 24, spatial distribution at 
3; = of electron normalized density pe - neKaQUcr) (line 1), pho- 
ton normalized density py = ny/(5aoncr) (line 2) and laser intensity 
pi = {E^ + B^)/E^ averaged over the laser wavelength (line 3) during 
the collision of two laser pulses with = 1-5 • 10^. 




Figure 8: At the instant of time t ^ 54, spatial distribution at 
_y = of electron normalized density pe = ne/(aoncr) (line 1), pho- 
ton normalized density py = ny/(5aoncr) (line 2) and laser intensity 
pi = {E^ + B^)/E^ averaged over the laser wavelength (line 3) during 
the collision of two laser pulses with <3o = 1.5 • 10^. 
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Figure 9: At the instant of time t ^ 62, spatial distribution at 
J = of electron normalized density pe = neliaQUcr) (line 1), pho- 
ton normalized density py = ny/(5aoncr) (line 2) and laser intensity 
pi = (E^ + B^)/Eq averaged over the laser wavelength (line 3) during 
the collision of two laser pulses with ao = 1.5 ■ 10^. 



y = 2000, a significant number of the electrons emit no 
photons. It follows from the obtained results presented 
in Figs. [Hand [21 that even in the limit ;^ <^ 1 when the 
regime of the photon emission is classical and the to- 
tal loss of the bunch energy due to the photon emission 
can be calculated in the classical approach the distri- 
bution function of the electron bunch after interaction 
cannot be described in the framework of classical elec- 
trodynamics if the mean photon number emitted by the 
bunch electron is small. 

We model the interaction between electron bunch and 
the laser pulse in quantum regime. The parameters of 
the laser pulse and the electron bunch are o-^ = 1.6, 
<2o = 100 and 70 = 2 x 10^. For these parameters;^ ^ 6 
and trad - 0.2. The electron distribution function after 
interaction with the laser pulse is shown in Fig.[3l It is 
seen from Fig. [T] that the bunch lost most of the initial 
energy. Moreover, some of the emitted photons decay 
and produce electron-positron pair. The similar eff'ects 
are observed in the experiments 1 15] and discussed, for 
example, in Ref. 1 16]. 

In this work we also present the results of 2D nu- 
merical simulations of QED cascade in the field of two 
colliding linearly polarized laser pulses. Laser pulses 
have the Gaussian envelopes at initial instant ^ = and 
propagate along x-axis towards each other. At initial in- 
stant of time the distance between laser pulses is equal 
to 2crx and coordinates of the pulse centers are x = -cr^ 
and X = (Tx. The cascade was initiated by 5 • 10^ 1 GeV 
photons situated near x = -cr^, y = 0. The values of 
simulation parameters are the following: A = 1.24 jim. 
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The spatial distributions of normalized electron and 
photon densities and laser intensity at three instants of 
time are shown at Figs.|4l[5]and[6l At initial stage of cas- 
cade development (Fig.|4l) e'e"" plasma density rapidly 
reaches the value of the relativistic critical density a^Ucr, 
where Ucr - m(j?l{^ne^) is the nonrelativistic critical 
density for e~e" plasma. The x-scale of created plasma 
is about the laser wavelength. Creation of dense plasma 
leads to substantial absorption of laser pulses and to 
the decrease of laser fields in the region occupied by 
the plasma. At later stages of the cascade development 
(Figs. [3 [51) e~e" plasma expands. 

Strong absorption of the laser pulses starts at t ^ 24 
when the plasma density becomes close to aoUcr- Due 
to the asymmetry in the initial position of the seed par- 
ticles one of the laser pulses partially has passed the re- 
gion occupied by the plasma ait ^ 24, while only small 
part of the front side of the other pulse has passed this 
region. Then the laser pulse that starts at x = cTx and 
the rear part of the other laser pulse are mostly absorbed 
by created e~e'^ plasma (see Fig.[5](c)). The asymmetry 
in the initial conditions also leads to asymmetry in the 
distribution of hard photons emitted by the electrons and 
the positrons. Normalized electron and photon densities 
and the distribution of laser intensity averaged over the 
laser wavelength at x-axis are shown on Figs.|7]l9l The 
asymmetry of the laser pulses absorption and hard pho- 
ton emission can be also seen in Figs.[7ll9l The steeping 
of the right edge of the electron distribution caused by 
momentum transmission from the absorbed laser pulse 
to e~e'^ plasma can be seen in Fig. [8] (line 1). 

In conclusion we present the simulation results 
obtained by the developed two-dimensional hybrid 
PIC/MC model. The model allows us to simulate QED 
processes in laser plasmas. First we analyze the QED 
processes during interaction of the relativistic electron 
beam with the intense laser pulse. It is shown that the 
QED eff'ects can be important even for moderately in- 
tense laser pulses when the number of emitted photons 
by single electron is not large. Then avalanche-like cas- 
cade development in the field of two colliding linearly 
polarized laser pulses is simulated. We show that over- 
dense (e~e'^) plasma is produced in the region where 
laser pulses overlap. A significant portion of the laser 
energy is absorbed as a result of production and heat- 
ing of the plasma. The anisotropic emission of hard 
photons is observed that can be used for development 
of bright radiation sources of gamma-quanta. The ob- 
tained results demonstrate that QED eff'ects can be ob- 
served in future laser facilities like ELI [ 20] and HiPER 
121], since the laser intensity (/ 2 • 10^"^ W/cm^) can 
be achievable in such laser systems. 
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